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Abstract 

We  develop  a  class  of  univariate  and  multivariate  life  distributions  for 
components  and  systems  which  operate  in  an  unknown  environment.  The 
environment  is  assumed  to  be  dynamic  and  is  described  by  a  shot-noise  process. 
Special  cases  result  in  a  univariate  distribution  with  monotone  hazard  rates,  and 
what  appears  to  be  a  new  family  of  multivariate  exponential  distributions  with  a 
singular  component. 
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1 ,  Introduction 


Consider  a  single  or  a  two-component  system  which  operates  in  an 
environment  that  may  or  may  not  be  identical  to  the  test  bench  environment.  The 
component  lifelengths,  when  assessed  under  the  test  bench  environment,  are  assumed 
to  have  failure  rates  X^Cu)  >  0,  i=l,  2,  for  u  ^  0.  We  shall  assume  that  the  nett 
effect  of  the  operating  environment  is  to  modulate  X-(u)  to  \j(u)  ?7(u),  i  =  l,  2,  where 
r?(u)  is  unknown  for  all  u  ^  0.  If  at  any  time  u,  the  operating  environment  is 
judged  harsher  (gentler)  than  the  test  bench  environment,  then  J?(u)  >(<)!; 
r?(u)  =  1  would  correspond  to  the  case  in  which  the  operating  environment  and  the 
test  bench  environment  are  identical.  Given  X|(u),  X2(u)  and  ^(u),  for  all  u  ^  0,  we 
shall  Judge  the  component  lifelengths  independent.  When  r?(u),  u  ^  0,  is  unknown, 
we  shall  describe  our  uncertainty  about  it  by  a  suitable  stochastic  process 
(ryCu),  u  ^  0}.  Under  such  circumstances,  there  is  an  induced  dependence  between 
the  lifelengths  of  the  two  components,  the  nature  of  which  depends  on  our  choice 
of  (mu),  u  ^  0).  The  case  X^Cu)  =  X^  i  =  l,  2,  with  r?(u)  ==  V,  and  uncertainty  about 
ri  described  by  a  gamma  distribution,  has  been  considered  by  Lindley  and 
Singpurwalla  (1986),  Currit  and  Singpurwalla  (1989)  and  the  references  therein;  the 
case  wherein  {r?(u),  u  ^  0}  is  described  by  a  gamma  process  has  been  considered  by 
Singpurwalla  and  Youngren  (1989).  In  this  paper,  we  shall  suppose  that  {r7(u),  u  ^  0} 
is  described  by  a  ‘‘shot-noise  process”,  a  motivation  for  which  is  given  below. 

Suppose  that  the  operating  environment  consists  of  a  series  of  events,  called 
“shots”  or  “jolts”,  whose  effect  is  to  induce  stresses  of  unknown  magnitudes  Xj., 
k=0,  1,  .  .  .  ,  on  the  components.  The  jolts  occur  over  time  according  to  a  Poisson 
process  with  a  known  rate  m(u),  u  ^  0,  and  suppose  that  whenever  a  stress  of 


magnitude  X  occurs  at  an  epoch  s,  then  its  contribution  to  r7(s+u)  is  Xh(s+u),  where 

the  att6nu3tlon  function  h(  )  is  positive  and  nonincreasing.  Thus,  if 

T  ^  T  <  .  .  .  ,  are  the  successive  epochs  at  which  stresses  of 

(0)  ^  (1)  (2)  ^ 

magnitudes  Xq,  Xp  Xy  ■  ■  ■  ,  respectively,  occur,  then 


r7(u)  =  ^  Xj.  h(u  - 
k-0 

with  h(u)  —  0,  for  u  <  0.  The  process  {5?(u),  u  ^  0}  is  called  a  shot-noise  process; 
see  Cox  and  Isham  (1980),  p.  135. 

As  an  example  of  the  above,  it  may  happen  that  an  unexpected  surge  of 
power  in  a  control  system  will  temporarily  increase  the  likelihood  of  failure,  but  the 
overload  often  decays,  lowering  the  likelihood  of  failure.  Another  example  pertains 
to  the  mortality  rate  of  individuals  suffering  a  heart  attack.  Since  the  heart 
muscle  repairs  itself,  the  rate  declines  with  elapsed  time  since  trauma,  but  the 
successive  bouts  of  attack  have  a  cumulative  effect  on  the  rate.  Situations  which 
necessitate  a  consideration  of  ‘‘cumulative  damage”  as  the  cause  of  failure  can  be 
incorporated  within  the  framework  of  shot-noise  processes  by  making  h(u)  a 
constant,  for  all  u  ^  0. 

When  X-(u),  the  Inherent  (or  base-line)  failure  rate  is  such  that  X^(u)  =  Xj, 
for  u  ^  0,  then  assuming  that  {n(u),  u  ^  0)  is  a  shot-noise  process  is  equivalent  to 
assuming  that  Xjr7(u),  the  failure  rate  of  the  i-th  component  under  the  operating 
environment,  is  also  a  shot-noise  process.  To  ensure  that  the  failure  rate  is  not 
annihilated  during  the  initial  phase  of  the  shot-noise  process,  we  require  that 


We  also  need  to  make  the  following  assumptions: 


Al.  Xj.  for  all  k,  where  “JI”  denotes  independence; 

A2.  X^JiXj,  for  all  with  Xj’s  having  a  common  distribution  G,  and 
A3.  m(t)  >  0  and  is  almost  everywhere  differentiable. 


2.  The  Survival  Functions 
2.1  The  Univariate  Case 

Let  L  denote  the  lifelength  of  a  single  component  system  whose  inherent 

failure  rate  is  X,  and  which  operates  in  an  environment  for  which  {r?(u),  u  ^  0)  is  a 

shot-noise  process  with  the  parameters  described  in  Section  1;  assume  that  Al  -  A3 

t  ,t 

hold.  Let  M(t)  =  m(u)du,  H(t)  =  h(u)du,  and  suppose  that  G*  is  the  Laplace 
transform  of  G;  then, 

Theorem  2.1.  Given  X,  M(t)  and  H(t),  F(t|X,  M(t),  H(t)),  the  survival  function  of  L 
at  T  ^0,  is 

.r 

G*{XH(t)}  •  exp  (-M{t)  +  G*{XH(t)}  m(T-u)du]. 

Lemoine  and  Wenocur  (1986)  outline  a  proof  of  the  above  theorem,  but  have 
overlooked  the  requirement  that  must  be  0.  The  proof  is  facilitated  by 

Lemma  2.2.  Consider  the  time  interval  [0,  r),  and  let  Tq^  <  Tq^  <  .  .  .  <  T^^^^  be 
the  epochs  of  events  in  a  nonhomogeneous  Poisson  process  with  rate  m(u),  u  ^  0. 
Then,  if  exists,  T-,  i  =  l,  .  .  •  ,  n  the  unordered  epochs,  are  independent  and 
identically  distributed  as  a  random  variable  U  whose  density  at  u,  o  u  <  t,  is 
m(u)/M(T). 

Proof  of  Lemma  2.2.  If  we  rescale  the  time  axis  by  M(t),  then,  on  the  new 
time  scale,  the  epochs  M(Tq))  <  M(T(2))  <  ....  are  Poisson  with  rate  1  [see  Cox 
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and  Isham  (1980),  p.  48).  Thus,  given  n,  MCT^^)  <  M(T(2))  <  •  <  are 

distributed  as  the  order  statistics  in  a  sample  of  size  n  from  a  uniform  distribution 
on  (0,  M(t)),  from  which  it  follows  that  the  unordered  epochs  M(Tp,  i=l,  .  .  .  ,  n, 
are  independent  and  uniform  on  [0,  M(t)).  Hence,  given  that  M  exists, 

P(M(Tj)  ^  u)  =  P(T-  M‘^(u))  =  u1M(t)  =  M  M"^(u)|M(t),  and  so 
P(Tj  ^  u)  =  M(u)1M(t),  from  which  the  result  follows. 

Proof  of  Theorem  2.1.  Suppose  that  n  jolts  occur  over  the  time  interval 
[0,  T)  at  epochs  T  =  (T^q)  <  <  ■••  <  with  each  jolt  inducing  a  stress  of 

magnitude  X.  Given  X,  n,  T  and  X,  the  survival  function  of  L  at  time  t,  is 

n 

exp  — |x  X  H(t)  +  X  X  H(t  —  , 

L  k=l 

which  because  of  A2  and  Lemma  2.2  can  also  be  written  as 

exp[-  {X  X  H(t)}]  [exp{-  X  X  H(t  -  u)}]". 

Unconditioning  on  X,  and  then  on  U,  the  above  becomes 

n 

G*(X  H(t))  I  G*fx  H(t  -  u))  m]^  du  . 

Jq  ^  "M(t) 

Our  final  step  is  to  uncondition  on  n  -  which  has  the  probability  mass 
function  (M(T))"/n!  -  and  to  simplify  the  ensuing  expressions. 
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2.2  The  Bivariate  Case 


Let  L^,  i  =  l,  2,  denote  the  lifelengths  of  a  two  component  system,  with  the 
i-th  component  having  as  its  inherent  failure  rate.  Suppose  that  the  system 
operates  in  the  environment  described  before.  Then, 

Theorem  2.3.  Given  Xp  X2,  M(t)  and  H(t),  F(Tj,  T2|X|^,  X2,  M(t),  H(t)),  the  joint 
survival  function  of  Lj  and  L21  for  0  ^  Tj  ^  T2j  is 


G*[XjH(Tj)  +  X2H(t2)]  •  exp  G*|^Xj  H(Tj— u^)  +  X2  H(t2— Uj)Jm(uj)duj 

L-'O 


exp  G*IX2  H(t2— U2)]ni(u2)du2  —  Mir 2) 


The  proof  of  this  theorem  is  facilitated  by  an  elaboration  of  Lemma  2.2. 
Lemma  2.4.  Consider  the  time  interval  (Tj,  T2),  and  let  <  •  •  < 

T/  \  be  the  epochs  of  events  in  a  nonhomogeneous  Poisson  process  with  rate  m(u), 

(n2' 

u  ^  0.  Then,  if  exists,  T^,  i=(nj+l),  ....  (n2),  the  unordered  epochs  are 
independent  and  identically  distributed  as  a  random  variable  U  whose  density  at  u, 
T,  <  u  <  Tn  is  m(u)/(M(Tn)  —  M(t,)). 

Proof  of  Lemma  2.4.  As  in  Lemma  2.2,  we  rescale  the  time  axis  by  M(t)  and 


observe  that  on  the  new  scale,  given  (n2  —  n^),  the  epochs  ^  ^^"^(nj+2)^ 

<  •••  <  M(T.  0  are  distributed  as  the  order  statistics  in  a  sample  of  size 
(n.,  —  nj)  from  a  uniform  distribution  on  [M(Tj),  M(t2));  the  proof  now  follows  via 
arguments  analogous  to  those  in  the  proof  of  Lemma  2.2.  II 


Proof  of  Theorem  2.3.  Parallels  that  of  Theorem  2.1,  with  the  additional 
proviso  that  for  a  Poisson  process,  the  number  of  events  in  [t p  T2)  is  independent 
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of  the  number  of  events  in  [0,  t^). 

Remark  2.5.  The  marginal  survival  functions  of  and  L2,  which  follow  from 
Theorem  2.3  by  setting  T2  and  Tj  equal  to  zero,  respectively,  agree  with  the  result 
of  Theorem  2.1. 

Remark  2.6.  A  generalization  of  Theorem  2.3  to  the  multivariate  (multi-component) 
case  is  straightforward,  with  the  k  >  2  variate  version  calling  for  (k-2)  additional 
parameters. 


3.  A  New  Bivariate  Distribution  with  Exponential  Marginals 


A  consideration  of  some  special  cases  of  Theorem  2.3  results  in  some 
interesting  distributions,  one  of  which  is  highlighted  here.  Suppose  that  for  all 
u  ^  0,  h(u)  =  1,  m(u)  =  m,  and  G  is  a  gamma  distribution  with  shape  a  and  scale  b. 
Then,  from  Theorem  2.3,  we  see  that  for  a  =  1,  the  joint  survival  function  of 
and  L2,  for  0  ^  T2  is 


-mb/(X,+X2) 

_  ^  (  b  1  f  1  ^ 

F(Ti,  T2IX1,  X2»  ==  U+XjTj+X2T2)^ 


\ _ b _ 1 

lb+X2(T2-Ti)J 


-mb/X 


2  -mr- 


(3.1) 


and  so  the  marginal  survival  function  of  L2,  for  T2  ^  0,  is 


mb 

Xo 

F  2(^2 1X2.  b,  m)  =  [  ^  j 


-  1 


-mr- 


(3.2) 
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Verify,  from  Theorem  3.1,  that  when  a  5^  1,  the  marginal  survival  function  of  L2, 


for  T2  ^  0,  is  of  the  form 


Ib+XrJ 


exp 


mb 


X(a-l) 


{> 


mr 


(3.3) 


If  we  assume  that  =  X2  =  X,  ®nd  that  m  —  X/b,  then  for  0  ^  ^  t-2> 

(3.1)  becomes 


FfTp  T2IX,  b,  m)  =  ^ 
and,  for  T2  ^  0,  (3.2)  becomes 


1 

-  mr^ 

+ 

mT2 

1 

+  mTj^ 

+ 

mT2 

(3.4) 


_  -mT2 

F.)(T2lX,  b,  m)  =  e 


(3.5) 


an  exponential  distribution.  The  results  for  T2  ^  and  F^tx^jX,  b,  m),  the 
survival  function  of  L^,  are  symmetric,  and  so  the  joint  survival  function  may  be 
written 


1  -  m  min(TpT2)  +  m  -  m  max(TpT2) 


F(Ti>  '^2!^’  ’3.  1  +  ni(Tj  +  T2) 


(3.6) 


Thus,  the  nonageing  property  of  a  component  in  the  test  environment  will  be 
preserved  in  a  shot-noise  environment,  if  the  stress  inducing  jolts  are  Poisson  with 
a  rate  X/b,  if  the  stresses  induce  a  cumulative  damage  on  the  component,  and  if  the 
magnitudes  of  stress  are  exponential  with  a  scale  parameter  b. 
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In  Figure  3.1,  we  describe  the  behavior  of  r2(T2),  the  failure  rate  function  of 
(3.2)  as  a  function  of  the  relationship  between  X2/'’  Verify  that 

r2(T2)  =  (X2/b)  (1 -t-mrl/d +(X2/b)T),  so  that  r2(T2)  m  as  T2  — »  o®,  and  that  r2(T2) 
is  increasing  (decreasing)  in  T2  and  asymptoting  to  m  depending  on  whether 
(X2/b)  <  (>)  m. 


Failure 

Rate 


(X2/b)  =  m 


Time  T2 


Figure  2.1.  Behavior  of  the  Failure  Rate  Function  of  F2(*) 


An  intuitive  explanation  of  the  above  behavior  of  the  failure  rate  function, 
as  a  function  of  the  relationship  between  (X/b)  and  m  is  given  in  Youngren  (1988). 
Of  particular  interest  is  the  reason  for  its  asymptoting  to  m.  Our  explanation  here 
is  that  under  a  shot-noise  process  with  cumulative  damage,  the  failure  of  a 
component  occurs  when  the  damage  exceeds  a  threshold.  With  the  inter-arrival 
times  between  the  jolts  being  exponential,  the  situation  here  is  analogous  to  the 
behavior  of  a  multicomponent  system,  with  component  lifelengths  being  exponential, 
for  which  the  life  distribution  is  a  gamma,  whose  failure  rate  asymptotes  to  a 
constant. 
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